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Maxwell's equations in an anisotropic inhomogeneous medium are transformed by means 
of the Stratton-Chu formula into a vector integral equation which couples the various elec- 
tric field components. In case the hypotheses of far-zone field and low frequency electro- 
magnetic waves apply, this vector integral equation can be approximated by a system of 
uncoupled scalar integral equations. This implies an approximate equivalence between the 
original vector integral equation and a system of modified scalar inhomogeneous Helmholtz 
equations. The conditions under which the system of uncoupled scalar integral equations 
can be solved by Neumann series are discussed, and the first three terms of the Neumann 
series are given explicitly. 



1. Introduction 

MaxwelPs equations for the electromagnetic field 
in an inhomogeneous anisotropic medium, such as 
the ionospheric magnetophisma, can be ^vl'itten in 
the form [Turner, 1954] 



curl E=— /xo 



V 

curl H=^^-= 



dH 



dD 

dt 



div D=0 



div H-0. 



(1) 



Assumhig harmonic time dependence e'"^ and that 
the ambient magnetic fiekl Bp is directed along the 
2-axis, one obtains the following vector wave equa- 
tion appropriate to an inhomogeneous anisotropic 
medium: 

curl curl E=coVoD=^2(v^ .g (2) 

y 

where A^o'^^Voeo and (k) is the conductivity tensor 
[Tiumer, 1954] of a magnetoplasma in a uniform 
magnetic field : 

/ B C 0\ 
Ck)=[-0 B \ (3) 

\ ^/ 

In the case of the ionosphere* A, B, and C are 
given by 



A=1-{X/U), B=l- 



-UX/(U'-Y'), 

C=iXYJiU'-Y') (4) 



where X, Y, and U=l—iZ are the usual definitions 
of magneto-ionic theory [Ratcliffe, 1959]. 



*For a cold magnetoplasma the collision frequency vanishes and definitions (4) 
reduce to 

A = l-X, B=l-X/(1-F2), C=iXF/(l-F2) (4') 



If the left-hand side of eq (2) is now expanded 
according to the usual vector identity, one obtains a 

term in grad div E=grad div [(k)"^ -D]. This term is 
a dyadic, and the various field components are thus 
inextricably coupled [Westfold, 1949] in an inhomo- 
geneous anisotropic medium. In all but the simplest 
cases the resulting systems of differential equations 
are practically intractable, and those formulations 
[Keller, 1953; Clemmow and Heading, 1954; Budden 
and Clemmow, 1957] which do exist are extremely 
involved, even in rectangular coordinates. It is the 
purpose of this paper to show how the problem may 
be simplified, at low frequencies and in the far-zone 
field, by using the Stratton-Chu formula [Stratton, 
1941; Silver, 1949] to transform to a vector integral 
equation ^ with uncoupled components. 

2. Transformation to a Vector Integral 
Equation 

In what follows, we shall assume that the medium 
is unbounded; i.e., consists of all of 3-space, so that 
there are no boundaries in the finite part of the 
medium. Silver-Miiller type radiation conditions 
[Silver, 1949; Mtiller, 1957] will be hnposed. The 
ionization will be assumed to be confined to a finite 
portion of space, so that at infinity A=l, B=l, and 
C=0. Thus at great distances, we have to do 
essentially with a free-space problem. 

The Stratton-Chu formula [Stratton, 1941; Silver, 
1949] reads as follows: 

r {G-VaVaE-E.VaVaG}(/F 

= [ {EaVaG-GaVaEJ-ik/Z. 

J bdryV 

(5) 



1 A somewhat similar formulation has recently been published by A. T. 
Villeneuve [1961] 
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In this expression n is the exterior normal to a region 
V of 3-space F3, and G and E are twice-diflerentiable 
vector functions of position over the closure of V, 
The vector E will of course be taken to be the total 
electric field, so that the left-hand side of the vector 
wave eq (2) may be introduced into (5). As men- 
tioned above, Silver-Mtiller type radiation condi- 
tions will be imposed upon E and H, so that the right 
hand side of (5) vanishes if the boundary of V is at 
infinity; i.e., y is all of 3-space. Following Stratton 
and Chu, we shall take G to be the free-space Green's 
function e'''^o^/R times a constant but otherwise 
arbitrary vector a. Thus 



For later use we note that 



(6) 



curl curl G=:a%+a47r5(P-Q) + (a-V)V^. (7) 

The result of substituting (2) and (7) into (5), 
and taking account of the radiation conditions, is the 
vector integral equation in the total field E : 

E(P)-^ J^ -l^klg{P,Q) (M).E(0 

+ {v'E)vg\dQ=0 (8) 



where 



(M) = (k)-(7) = (<.)/(zW„), 



(9) 



and (I) is the unit dyadic. Since the operator V 
applies to the source point Q, we may also write the 
integral equation (8) as 

E(P)=^J^ g{P,Q) ■l^kl{M).E(Q) 

+(ih+j^^k(v'E}ydQ, (10) 

A 

where R denotes a unit vector in the Yp—Yq direction. 
If we write the total field as the sum of incident and 
scattered fields : 



E(P)=E^-(^)+E^(^), 



(11) 



then the integral equation (10) may be written as 
the inhomogeneous integral equation in the unknown 
scattered field : 



E^ 



(P) = r(P)+^ J^^^(P,e) ^kl{M)^E^{Q) 

+(^^o+;|)R(V-Eo]-rfe, (12) 

where F'(P) is a known vector function of position 
defined in terms of the incident field by the formula : 

+(^^0+1) R(V.Eo|rfQ-E^-(P). (13) 



3. The Low-Frequency, Far-Zone Approxi- 
mate Form of the Integral Equation 

In the sequel attention will be restricted to the 
case of low frequency electromagnetic waves, i.e., 
those fields for which the free space wave number is 
such that 

^o«l. (14) 

However, this hypothesis does not imply that the 
first term of the integrand in (10) or (12)) is neces- 
sarily small, for even in a weakly ionized magneto- 
ionic medium the product ||A:g(if)|!==coAto||(o-)Tl can 
be large whenever co<^<Caj^<^<^coe/w. 

We further require that the field point P be in 
the far-zone with respect to the distribution of 
ionization, the latter being assumed to be confined 
to a finite volume Vc of space. Thus we can impose 
the hypothesis that 



Pi>Q»^o-=Xo/(2x) when QeV,. 



(15) 



The second term in the integrand of (10) (and 
(12)) contains the quantity div E, which since we 
assume quasi-electrostatic equilibrium, can be iden- 
tified as proportional to the polarization charge 
density : 

V.E = pp/eo 

But the equation of continuity implies, since 



J=^— ; that 



V- E= 



V J 



icoeo 



(16) 



where J is a conduction current occurring only in 
the ionized portion of space : 



That is 



V- E 



^{0: 



i={<j) .E. 
V((M).E)inFe 



in V,-Vc. 



(17) 

(18) 



For those source points Q in V^ the presence of 
ionization implies that (o')p^(O), and so (M)f^(O) 
and V-Ef^O; for all other Q, (o-)=-(0), and so (M) 
and V-E vanish. Therefore the only contribution to 
the integrals in (10) (and (12)) can arise from those 
Q which are both in Vc and such that P>>Xo. If 
we now impose the hypothesis that 



1V((M).E)|<|(M).E|, 



(19) 



(i.e., a slowly- varying distribution of ionization), 
then (19) together with (14) and (15) implies that 
the second term in the integrand of (10) (and (12)) 
is of higher order than the first. We may therefore 
approximate (10) by the polar integral equation 



^iP)-%^yjiP,Q)W-^iQ)dQ- 



(20) 
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We state this result as a theorem : 

Theorem 1. Under hypotheses (14), (15), and (19) 
the vector integral equation reduces to the approxi- 
mate form 

E(P)=§J^^ g{P, Q)iM).EiQ)dQ. (20) 

Corollary 1.1. Under the same hypotheses the in- 
tegral equation (13) in the scattered field E^ reduces 
to the approximate form 



E 



'(P)=P(P)+§: J^ g{P,Q) (M) . WiQ) dQ, 



(21) 



where F*(F) is given by (13). 

Corollary 1.2. For a cold plasma {v=Q) hypoth- 
esis (19) may be replaced b3' the hypothesis that 
X>1, F»l over Vc. 

Proof. Since 



||M|| = 



X . XY 



i^Y'2 "i-r^ 







—I 



. XY 



X 



, 



l-Y' l-Y 
-XJ 

=x(2,^ + l)=A'(^ + l) (22) 

the conclusion follows immediately. 

4. Transformation to a System of Uncoupled 
Scalar Integral Equations 

In an appropriate orthogonal basis (ei, 62, 63) 
(which may, in particular, be the base vectors for 
cartesian or cylindrical coordinates), 

(.l/).E = [(B--l)^i+C'^2]ei 

+ [{B-l)E,-CE,]^2^-{A-l)E,e,. (23) 

The vector integral eq(20) can then be written as 
the coupled system of scalar integral equations 

E,{P)^^^j^^g{P,Q)[-CEm + {B-l)E,]dQ 
^3(P)=2 J^^ g{P,Q){A-\)EmdQ 



(24) 



with <a corresponding form for eq(21). 



Now define the following quantities: 

cfl,2(P)=£'l±'i£'2 

<^z{P)=E, 

^.2iP)=Fi±iF2 

MAP)=B-lTiC 

j/iP)=A-l 



(25) 
(26) 
(27) 
(28) 
(29) 



We note that Si and S2 represent counter-rotating 
circularly polarized components about the Xa-axis, 
while (03 is the longitudinal field component directed 
along the a:;3-axis. Further, under the hypothesis 
(X>1, Y>1), of corollary 2 



and 






s/=-X=-<4/^''. 



(30) 



(31) 



Then the system (24) can be written as the follow- 
ing uncoupled system of scalar integral equations 

in X'' /^ A^ ' 
ni (b 1, & 2y <OZ' 



V. 47rjt3 



(32) 



with a corresponding uncoupled fonn for the com- 
ponent form of eq(21) : 

(33) 

5. Equivalences Between the Vector Wave 
Equation and Approximate Systems of 
Uncoupled Scalar Helmholtz Equations 

The following equivalences can now be stated: 
Theorem 2. Under hypotheses (14), (15), and (19) 
above, the vector wave eq(2) can be approximated by 
the uncoupled system of scalar Hehnholtz equations : 



/Ei+iEA /B-iC 








AE,-iEA^kl[ 


Bi-iC 





\ E, J \ 





A 



(E,+iE\ 

Ei-iE^ =0. (34) 
^ E, J 
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Prooj. This follows immediately from applying 
to the system (32) the equivalence between the 

integral equation u{P) — K g(P,Q)p(Q)u{Q)dQ=0 

and the boundary value problem consisting of 
V^u+A[l + p(P)] u{P) = together with homoge- 
neous Dirichlet or Neumann boundary conditions 
and/or a radiation condition, and at the same time 
taking into account definitions (25), (26), (28), 
and (29). 

Theorem 3. Under hypotheses (14), (15), and (19) 
above, and the assignment of an E*(P) such that 
there exists an/^^(P) satisfying 



=r'(P)=^J^ giP,Q)/WdQ, 



(35) 



the vector wave equation curl curl E^—kl(K)-'E^ 
=P(P) is essentially equivalent to the uncoupled 
S37^stem of scalar inhomogeneous Helmholtz equations : 



/El+iEl\ /B-iC 








vHEi-iEi]+kU 


B+iC 





\ El J \ 





A 



/Et+iEi\ /ji+ifh 

lEl-iEl\=Ut-ifil (36) 

Proof. The theorem follows in the same way 
as theorem 2 from the equivalence between the 
integral equation 

u(P)-AJ9{P,Q)p(QMQ)dQ=jg(P,Q)q(Q)dQ 



and the inhomogeneous Helmholtz equation 
V^u-{-A[l-{-p(P)] u(P) = Cj[{P) with homogeneous 
boundary conditions and/or a radiation condition. 

6. Solution of the Uncoupled Integral 
Equations by means of Neumann Series 

Each of eqs (33) is of the form of a singular 
Fredholm integral equation with symmetrizable 
kernel. Since the domain V^ has been taken to 
be all of 3-space, and so is infinite, the integral 
equations will be singular. However, Neumann 
series solutions of the integral equations can still 
be defined [Schmeidler, 1955] in the usual form and 
subsequently shown to actually converge to the 
solutions of the respective integral equations, pro- 
vided the kernels Ki,2{P,Q)=g{P,Q).^u2(Q) and 
K^{P,Q)=g(P,Q)s^(Q) are integrable over F3 with 
respect to every function d^ (P) in L2 and with 
respect to every function drjiQ) in L2, and both of 
the following iterated integrals exist and are equal : 



j^^'{P)(^j^Kn(P,QW{Q)dQ^dP 

The way to proceed is described by the following 
theorem [Schmeidler, 1955]: For a*^kernel K(P,Q) 
which satisfies the three requirements just stated 
and whose associated bilinear form is absolutely 
bounded: 

I r ^^ KiP,Q)^'*{PW(Q)dPdQ\<M, 

The Fredholm equation of the second kind 

y(P)~AJ^K{P,Q)y{Q)dQ=F(P), yeL, 

has a unique solution in the form of a mean-square 
convergent Neumann series 

yiP)=l.i.m.-^F{P)+AJKiP,Q)F(Q)dQ 

+A'JK(PQ) (^JK(Q,Q')F{Q')dQ'yQ+ ■ ■ X- 

In the present application the hypotheses of this 
theorem are equivalent to 



I r „ giP,Q)mA(M'*iPW{Q)dPdQ 



<M,,2 



and 



(37) 



iJ. 



„ ^„. g{P,Q)j./{Q)r{PW{Q)dPdQ 



<M^ 



(38) 



where ?'(P) and ri'{P) are arbitrary L2 functions on 
Vz, and 



M,<irK (fif= 1,2,3) 



(39) 



The zero-order terms in the Neumann series solu- 
tion for S\{P), (fl(P), and (fl(P) are simply ^l(P), 
.j^iiP), and Fl(P), respectively. The respective 
first-order terms (Born approximation) are given by 
the integrals : 



A.0 

47r. 



and 



f g{P,Q)^^UQ)^nAQ)dQ (40) 

^^j{P,Q)s^{Q)Fm)dQ. (41) 
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The third terms in the corresponding Neumann 
series are, respectively, 



and 



(S)X '^^^^^' ^)^'-^^^X^ ^^^' Q')-^.,2(Q').n2(«'MQ' 



(42) 



(43) 



etc. Since A:o/47r«l in the low frequency case, the 
first two or three terms in the Neumann series will 
ordinarily suffice for numerical purposes. 
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